train_set=fom1.in
order=4
fom_model=fom.model
vom_model=vom.model
likelihood=fom1.ll
eta=0.2
simulate=vom.simul.in
n=1000
length=90
fom_simul_model=fom.simul.model
vom_simul_model=vom.simul.model
likelihood_simul=fom1.simul.ll
./vom --fom --input=$train_set --order=$order --output=$fom_model
./vom --grow --input=$train_set --model=$fom_model --output=$vom_model --eta=$eta --prune
./vom --simulate --model=$vom_model --output=$simulate --n=$n --length=$length

./vom --fom --input=$simulate --order=$order --output=$fom_simul_model
./vom --grow --input=$simulate --model=$fom_simul_model --output=$vom_simul_model --eta=$eta --prune

./vom --likelihood --input=$simulate --model=$vom_model --output=$likelihood
./vom --likelihood --input=$simulate --model=$vom_simul_model --output=$likelihood_simul

diff $likelihood $likelihood_simul | less